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ABSTRACT 

We present an investigation of the multiwavelength data on HESS J 1303-631, an unidentified TeV source 
serendipitously discovered in the Galactic plane by the HESS collaboration. Our results strongly suggest 
the identification of this particular source as the remnant of a Gamma-Ray Burst (GRB) that happened some 
few tens of thousands years ago in our Galaxy at a distance on the order of > lOkpc from us. We show 
through detailed calculations of particle diffusion, interaction and radiation processes of relativistic particles 
in the interstellar medium, that it is possible for a GRB remnant (GRBR) to be a strong TeV emitter with 
no observable synchrotron emission. We predict spectral and spatial signatures that would unambiguously 
distinguish GRBRs from ordinary supernova remnants, including: (1) large energy budgets inferred from their 
TeV emission, but at the same time (2) suppressed fluxes in the radio through GeV wavebands; (3) extended 
center-filled emission with an energy-dependent spatial profile; and (4) a possible elongation in the direction 
of the past pair of GRB jets. While GRBRs can best be detected by ground-based gamma-ray detectors, the 
future GLAST mission will play a crucial role in confirming the predicted low level of GeV emission. 
Subject headings: cosmic rays — diffusion — gamma-rays: bursts supernova remnants 



1 . INTRODUCTION AND OUTLINE 

So far, GRBs have only been identified at cosmological dis- 
tances. It is believed that GRBs are caused by highly relativis- 
tic outflows with bulk Lorentz factors V > 100 that form a 
pair of opposite jets. Even after correcting for narrow beam- 
ing, the energy radiated by long-duration (> 2 s) GRBs typ- 
ically is ~ 10 51 erg, and the estimated kinetic e nergy of the 
jets reaches values of ~ 1 52 erg (jMeszaros 2002; Bloom et alJ 
120031 iBereer et al.ll20 04). Relativistic shocks convert (accel- 
erate) all particles from the incoming plasma they encounter 
to relativistic energies by randomizing their velocities in the 
comoving frame. In the sta tionary frame the mean ener gy per 
particle is E ~ m p c 2 T 2 /2 (Blandf ord & McKed [l976l). The 
total energy of cosmic rays (CRs) produced by a GRB can 
thus be up to ~ 100 times higher than the ~ 10 50 erg pro- 
duced by no nrelativistic shocks of typical supernova rem- 
nants (SNR) ( Berezinskv et al. 1990). GRBs are prime can- 
didate sources for the extragalactic component of CRs, and 
they m ay also be sou rces of ultrahigh energy CRs in our 
Galaxy (iDerme Bum. Calculations taking into account the 
star formation history and the GRB beaming factor predict 
a Galactic rate of one long GRB every 30 00 to 10 s years 
JWicket alJkOOtt iDermer & Holmes! 12005). A GRB that 
occurred ~ 1 Myr ago at < 1 kpc from us could expl ain the 
knee in the observed CR spectrum at E > 10 15 eV dWick et al.l 
2004). Ohers have suggested that giant radi o shell objects 
in the Milky Way, such as the SNR W49B Jsee Ioka et al.1 
2004) or kpc -size HI super-shells (Efrem ov et alJl9 98) could 
be remnants of GRBs. 

The HESS (High Energy Stereoscopic System) collabora- 
tion has recently dis covered a population of 7-ray sources 
in the Galactic plane (lAharonian et al.ll2005bl 1200 6) that are 
TeV-bright but are often quiet in all other wavelengths and re- 
main unidentified. All previously detected TeV sources show 
strong X-ray emission which is the synchrotron counterpart of 
the Inverse Compton (IC) TeV radiation from multi-TeV elec- 
trons. The suppression of low-frequency flux is unexpected 



within conventional models. Mitra (2005) has proposed that 
this feature points to radiation production in the deep grav- 
itational potential wells of black holes formed by Galactic 
GRBs. However, the extended appearance of these sources 
seems to contradict this explanation. 

The estimated rate of Galactic GRBs implies a likelihood 
of between one and several "young" (age t > 10 4 yr), GRB 
remnants in the interstellar medium (ISM) at distances up to 
af~ 10-20 kpc from us. Some > 10 4 years after their acceler- 
ation by the relativistic shocks of GRBs, multi-TeV particles 
would diffuse over distances of about > lOOpc. Therefore, 
if detected in TeV 7-rays, Galactic GRBRs should appear as 
large plerionic ('center- filled') nebulae in the ISM. 

While probably most of the unidentified TeV sources would 
eventually be associated with conventional Galactic objects, 
here we argue that at least one source, HESS J 1303-631 
( Aharonian et al . 2005a), shows features that strongly suggest 
its identification with the remnant of a GRB that happened 
> 10 4 (-k//c) years ago in our Galaxy at a distance d~ (10- 
15) kpc from us. We show that this can explain both the total 
CR energy in the source, and its non-detection below TeV en- 
ergies. Furthermore, we discuss observational signatures that 
can be used to identify GRBRs unambiguously. 

2. INVESTIGATION OF THE TEV ENERGY SPECTRUM 

Among the unidentified TeV sources discovered by HESS, 
J 1303-631 is the brighte st source with the most detailed data 
l lAharonian et al.l2005al) . The measured (0.38-12) TeV differ- 
ential spectrum, approximated as a power-law F(E) oc E~ a i , 
has a mean spectral index a 7 « 2.44 and an integrated en- 
ergy flux /e ~ 1.87x 10 _11 ergcm" 2 s _1 . So far, no emission 
has been detect ed in any other wave ba nd, including the X- 
ray band (iMukheriee & Halpernl 12005). The source is ex- 
tended and is approximated, within current statistical errors, 
by a spherically symmetric 2-dimen sional Gaussi an distribu- 
tion with inherent width 9 src = 9'. 6 (Aharonian et al. 2005a). 
The center-filled appearance of the source and the gradual 
steepening of the TeV spectrum suggest the emission origin 
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from a population of high-energy particles injected into the 
ISM from a point source and evolved by energy-dependent 
diffusion. 

Absence of a detectable synchrotron flux implicitly sug- 
gests a hadronic origin of the TeV radiation. Since the two 
7-rays from a 7r°-decay each carry a fraction of the inci- 
dent proton energy of about E 1 ~ 0.075 E p , the observed 
flux requires U p ~ And 2 fEt pp r]Q l total energy in protons with 
E p > 5 TeV. Here, t pp = (K pp a pp cn H )~ l is the /?/;>-cooling time, 
and rjo — 1 /3 is the fraction of the lost proton energy that is 
converted into 7r°-mesons. Substitution of the cross-section 
(Tpp ~ 40mb and the inelasticity K pp ~ 0.5 gives 

U p ~\.2y.\0 A9 dl vc n^(\+S)exg. (1) 

Here dk pc =d/ 1 kpc, and S accounts for protons with E < 5 TeV. 
For standard SNR acceleration, a broken power-law spectrum 
with a 2 ~2.44 above and a,\ =2.0 below 5 TeV gives the 
smallest possible 5«4.9. A more plausible SNR spectrum 
with c\!i =2.2 predicts S~2A. 

Integration of the proton energy distribution 
forme d by d iffusion from a point 'impulsive' source 
( Atov arTet al.l 119951) along the line of sight gives 
n(E, 0,f) = Nq{E)-k-H^ { exp[-C(£, 0)/2], where l m = ^2D(E)t 
is the mean diffusion length (along each of the axes) 
for a particle with energy E, diffusion coefficient 
D(E) = D 27 1 27 (£'/10GeV)*cm 2 s- 1 in a power-law ap- 
proximation, and ((E,9) = (Od/ldid 2 - Integrating over 
a circle of angular radius S , we derive a proton energy 
spectrum N(E;9 < 9 S ) = N (E) [1 -exp(-C(£\ f )/2)]. 

The observed spectral index a 1 = 2.44 is significantly 
steeper than at\ « 2. 1-2.2 expected for a typical SNR. We ex- 
plain this steepening as an unavoidable result of energy de- 
pendent diffusion. Given the HESS point spread function 
(PSF), the derived N(E ;9<9 S ) suggests Co = C( 1 Te V, 9'.6) ~ 
1. Calculations nicely fit the data for Co = 0.61, and ex- 
clude Co > 0.8 or Co < 0.45 (see Fig. 1). This defines the fol- 
lowing key-relation between the source distance and its age 
t = ?kyr 10 3 yr (we count t from the time when the first light 
from the source reached us, and neglect the source expansion 
during the light crossing time 2/dif/c <C t) : 

d kpc = 0.935 x 10 3S /\D 27 f kyr ) 1/2 Co /2 • (2) 

This constrains the energy required to explain both the spec- 
tral and the spatial appearance of the source: 

U p « 1.1 x 10 49+35 (1 +S)D 27 f kyr «H 1 Co erg. (3) 

The diffusion coefficient in the ISM should have 5 ~ 0.5 to 
explain the CR spectral index qck = ot\ + 5 sa 2.7. A typi- 
cal value of D(E) at lOGeV is D 21 > 10 jBerezinskv et all 
119901 iPtuskin & Soutoull 119981) . A lower limit of D 21 > 
0.3 can be derived from the grammage X = m p n H cto — 
10 ge m" 2 that explains the ob served composition of CR nu- 
clei (le. g.. Buckley et al.l ll994). and does not allow CRs to 
stay for more than to ~ 10 7 yr in the Galacti c disk at h eights 
%sk < 150pc where n H > 0.5 cm" 3 jPtuskin & Soutoull 1998b . 
Values of D 2 j much smaller than the typical ones for the ISM 
would imply an enhanced level of turbulence and can be ex- 
cluded, as the associated enhanced magnetic field would re- 
sult in radio emission exceeding the flux upper limit derived 
below. 

With these results, we are now in the position to argue 
that the source is a GRBR and not a typical SNR. Using 
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FIG. 1. — Observed (stars) and modelled (lines) angular distributions of 
7-rays from HESS J1301-631 (normalized at 8 = 0). From bottom to top, 
the lines show the distributions for Co = 1-22, 0.862, 0.61, 0.43, and 0.305, 
convolved with the PSF oc aexp|-0.5(fl/2o-i ) 2 1 +fcexp[-0. 5W2g] ) 2 1, with 
parameters u\, ai and a/b derived from the PSF plot in Aharo nian et alj 
1 2005a). The full dots show the best-fit to the experimental data in the form 
of PSF-convolved 2D Gaussian distribution with width # slc = 9.6'. 



D 2 i ^ 0.3, we constrain the source age to ?k yl - < 0.2t/^ pc . 
Even for a SNR with a\ =2.0, an age t k yr < 1.8, and loca- 
tion at c/^3kpc in the ISM with «h ~ 1cm" 3 , equation (1) 
predicts U p > 6.4 x 10 50 nH erg. For a typical 10% CR in- 
jection efficiency by a SNR, this requires a kinetic energy 
t4in ^6x 10 51 erg. In principle, such energies cannot be ex- 
cluded for some extremely powerful supernovae (SNe) of type 
Ib/c or Iln (which in fact might be GRB s) that constitute 
< 1% of all SNe k.s. Sveshnikoval 120031) . This implies a 
total of A^sn ;$ (0-l-0.3)fk yr such events in our Galaxy. The 
expected number of such SNRs at a distance d kpc « 3 is only 
£~ jV SN (<4 pc /15) 2 < (0.7-2.1) x 10" 2 . Furthermore, the non- 
detection of the shell of such powerful, close and young SNR 
in the radio and X-ray bands is difficult to explain. 

An assumption that the source is an SNR insid e a molec- 
ular cl oud with /in = 56 cm" 3 at d = 2.1 kpc ( Aharoni an et al.l 
2005a) relaxes the extremely high SNR energy requirement 
and increases £ to acceptable values. However, the problem of 
the missing shell of this still powerful and young SNR in such 
a dense medium persists. Even for a total energy in relativistic 
electrons of about 1% of the protons, which is a typical lower 
limit for an SNR, the synchrotron flux would significantly ex- 
ceed the flux upper limit at 5 GHz (see Fig. 3). 

For a source age of t > 10 4 yr, equation (|2ji moves the source 
to d > 10 kpc and increases dramatically the energy require- 
ments for an SNR that accelerates particles by non-relativistic 
shocks. The full radius of the source (detected up to 9 <25') 
increases to >70pc implying that it is located either in or- 
dinary ISM or entirely inside (to have a quasi-uniform profile 
for «h) a giant molecular cloud. Figure 2 shows the fluxes cal- 
culated for «h= 1 cm" 3 , d= 12kpc and t = 1.5x 10 4 yr. Two dif- 
ferent proton energy spectra are considered. One implies an 
'extreme SNR', with a2=2.2 at E>E\,± =5 TeV but a\ =2.0 at 
E <£brk- The total energy derived is L^.snr= 8. 1 x lO 51 ??^ 1 erg. 
The 'GRBR' spectrum assumes a\ = 1 .5 implying protons ac- 
celerated by a relativistic shock. For qualitative estimates the 
break energy could be defined as £brk^0.5r 2 GeV, where Y\ 
is the shock Lorentz-factor at the time when most of the initial 
shock energy is transferred to relativistic protons (see Sec. 3). 
Fitting the data we find T\ < 100. Thus T\ can be as high 
as the initial T > 100 of GRB jets JMeszarosll2002l) . Figure 
2 assumes Ty = 100. This reduces the total energy required 
for the GRBR model to C/ p ,grbr = 4.1 x 10 51 «h erg. Values 
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FIG. 2. — Hadronic (7r°-decay) model of the 7-ray emission from HESS 
J1301-631, calculated for a GRBR of age t = 1.5 X 10 4 yr at d = 12kpc in 
the ISM with «h = 1 cirT 3 . Calculations assume a broken, "GRBR type", 
power-law proton spectrum, with cq = 1.5 below £brk = 5TeV. The thick 
dashed line shows the predicted energy flux when all emission within one 
degree from the source center is considered. The spectrum is considerably 
harder than the measured one (data points and the bar). Owing to technical 
difficulties associated with the determination of the energy spectrum of an 
extended source, the H.E.S.S. collaboration derived the spectrum using only 
photons detected within the central 13.4' radius and scaled it according to 
the total observed flux. We predict that the energy spectrum derived using all 
the TeV photons will follow the dashed line. The solid line shows the flux 
from the intrinsic 8 < 10' region (this effectively translates to 13.4' size af- 
ter correcting for the PSF). The dashed-dotted line and the dotted lines show 
the fluxes from 10' -25' and a 25'-l° annuli, respectively. The 3-dot-dashed 
curve is the flux from the inner 10' angular region for an SNR-type spec- 
trum of CRs, a, = 2.0. Upp er flux limit of EGRET for HESS J1303-631 
1 Mukheriee & Halpern 2005) and the 1-yr flux sensitivity of GLAST are also 
shown. 



of £t,rk < 5 TeV, which are not excluded, would lead to an 
increase U ^.grbr oc (STeV/Ebrk) 0,2 (for £brk > 0.1 TeV), and 
even slower for U p .snr- 

The HESS collaboration dis cussed the giant (r ~ 1 44pc) 
molecular cloud G303. 9-0.4 (iGrabelskv et alJ [1988) with 
«h — 10-20 cm" 3 located at 12 kpc from us as a possible 
host of the source. For the SNR-type spectrum, the mini- 
mum CR energy U p .snr =(4-8)xl0 50 erg is problematic even 
for a very powerful SNR. Furthermore, sub-relativistic pro- 
tons alone would contain sufficient energy for heating (at a 
rate ~3 x 10 35 ergs _I ), and possibly also disintegration, of the 
cloud via Coulomb losses. Because of the "low-energy cut- 
off" expected below E\,±, these constraints do not apply if 
HESS J 1303-631 is a GRBR. It is also possible that the cloud 
is merely a line of sight coincidence. Even at a distance of 
d =15 kpc a GRBR would still be able to power the inferred 
proton energy density in the ISM. 

Leptonic models require less power than hadronic models. 
IC radiation of multi-TeV electrons on the CMB and Galactic 
far-infrared photons could reproduce the observed spectrum. 
A broken power-law electron spectrum with a e .\ =2.1 below 
10 TeV, and a e .2=3.2 above leads to the minimum required 
energy U e ~2.4x 10 48 afj; pc erg. However, the synchrotron ra- 
dio flux of the leptonic model exceeds the flux upper limit that 
we deri ved from the Parke s -MIT-NR AO southern 4850 MHz 
survey (Condon et al] ll993l) even for an ISM magnetic field 
fiiSM = 3/iG (see Figure 3). We estimate the limit by measur- 
ing a mean flux of 0.0032 Jy/beam in the 9'. 6 radius region 
centered on the source. Calculating the flux in background re- 
gions of the same aperture, we estimate a standard deviation 
of 0.0033 Jy/beam. Using a beam width of 7' FWHM, we 
then obtain a 5a upper limit of 0. 12 Jy on the flux from within 
9'. 6 from the source centroid. 

We have also computed the IC fluxes from the sec- 
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FIG. 3. — Synchrotron fluxes (shown in the main frame) expected from 
HESS J1301-631 in the case of a leptonic (IC) origin of TeV fluxes (shown in 
the inset). The fluxes expected from angular distances 9 < 10' and 9 > 10' 
are plotted with heavy solid and dot-dashed lines, respectively. The total IC 
flux (dashed line) is normalized to the observed flux at 1 TeV. The source age 
t= 10 4 yrandB = 3A*G are assumed. The thin solid curve shows the syn- 
chrotron flux from accelerated (primary) electrons with total energy 1% of 
relativistic protons, to be expected from a SNR at d = 2. 1 kpc in dense cloud 
with «h = 56cm~ 3 i Aharonian et al. 2005a) and enhanced level of turbulence 
and magnetic field, B = 10 /iG, where the observed TeV flux would be ex- 
plained by the protons. The synchrotron radiation from secondary electrons 
is shown by the full-dotted line. 



ondar y ir -decay e lectrons and from /3-decay electrons 
(see Ioka et al. 2004), accounting also for their radiative en- 
ergy losses. Both fluxes are almost two orders of magnitude 
lower than the 7r ( )-decay flux. The ratio of synchrotron to IC 
fluxes depends only on the ratio of the magnetic field to the 
target photon energy densities. Hence the radio flux from the 
secondary electrons for the GRBR model in Figure 2 does not 
violate the derived upper limit in so far as their IC emission is 
well below the observed TeV flux. 

3. DISCUSSION AND PREDICTIONS 

A GRBR should differ from a SNR not only by the energet- 
ics, but also by three characteristics of its relativistic particle 
population. First, the prompt differential spectrum of protons 
(and electrons) behind the relativistic shock cannot be a single 
power-law N(E) oc E~ a with a > 2.0 as in the case of non- 
relativistic shocks. As shocks with Lorentz-factor T{t) 3> 1 
boost all particles they encounter to relativistic energies, the 
prompt spectrum should break to an index ai iPr0 mpt < 1 below 
some energy E^ A <m p c 2 T 2 (t)/2 (observer frame) to satisfy 
the laws o f energy a nd particle number conservation across 
the shock (Katzl [19941) . The final power-law index ckj (pre- 
sumably < 1 .5) and the break energy E\,± are formed while the 
shock decelerates, and will depend on a number of processes, 
such as diffusive escape of relativistic protons from the shell, 
their adiabatic cooling and reacceleration, etc. At E > E\,± a 
spectra l index a? ss 2.2 is ex pected to form in Fermi-type pro- 
cesses (Kirk & Duffvi fl99^) . For the current data, £t,rk could 
be as high as 5 TeV (in Figure 2), but it could also be smaller, 
leading to GeV 7-ray fluxes detectable by GLAST. 

The second feature of GRBRs that we predict is that a 
substantial fraction of the initial shock energy is eventually 
transferred to high-energy protons, while relativistic electrons 
would contain only a tiny fraction of it at times of non- 
relativistic evolution of the GRBR shell. A GRB shock be- 
comes non-relativistic, and als o quasi-spherical, o n timescales 
< lyr at distances < 0.1 pc ( Livio & Waxman 2000k At 
this stage, magnetic fields of ~0.1G ( Berge r et alJl20041) are 
still high and will (synchrotron) cool electrons to a few GeV 
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within ~lyr. As the initial kinetic energy (5-10)xl0 51 erg 
of the jets is much larger than the radiated (and the mag- 
netic field's) energy, by then most of it should be in rela- 
tivistic protons. The evolution of a non-relativistic shell that 
is dominated by the relativistic energy, and not by the ki- 
netic energy of the shock, requires a fully relativistic treat- 
ment and remains unexplored. Diffusion of relativistic pro- 
tons (thus escaping adiabatic losses) upstream may broaden 
the shock precursor. Because of the preferential accelera- 
tion of_rwticJes_wjfli_I i orejitz_factors larger than the precur- 
sor (Berezhko and Ellison 1999), energetic protons will read- 
ily be accelerated while acceleration of electrons will be sup- 
pressed. This would also apply to the spherical "supernova"- 
type component (which, however, may not be developed by 
GRBs at all, see MacFa dTen & Wooslevl (119991) 1 as this non- 
relativistic shock has to pass through the central < lOpc re- 
gion dominated by the energy of relativistic protons from the 
(quasi-spherically opened up) GRB jets. We also note that 
in the presence of relativistic protons, nonrelativistic shocks 
may evolve (and dissipate) much faster than predicted by the 
standard Sedov-Taylor solution for SNR shocks propagating 
in "cold" plasmas (Toptvain 2000). 

Third, considerable energy, ^(10-20)%, of the acceler- 
ated protons is c hanneled i nto collimated beams of neutrons 
jAtovan and Dermerl2003l iDermer & Atovanl2003l) from in- 
teractions p+7^«+7r + with X-ray photons from the GRB. 
First, this mechanism assures a minimum level of total energy 
of the protons that avoid adiabatic losses in the shell, U„ > 
10 51 erg. Furthermore, neutrons with energies E n will /3-decay 
into protons on length scales of l n (E n ) ~ (£'„/100TeV)pc. A 
characteristic signature of a GRBR might then be an elonga- 
tion of the TeV 7-ray emission due to conversion of the ini- 
tial double-sided neutron beam into protons, which will sub- 
sequently be injected into the diffusion process along the jet 
path on length scales ~ l-10pc (for£„ < IPeV) . Depending 
on the contribution of neutrons to the high-energy proton pop- 
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ulation, the two-dimensional Gaussian profile may have an as- 
pect ratio of up to A = (an — o"j_)/(<7||+c_l) with a± = ldif{E,t) 

and <j|| = yj l^(E) + l^ i{ (E ,t). For a source size Zdif }t 10pc this 

would predictA(£')^(0. 01-0.1). However, neutron beams can 
also affect the shape of the GRBR nebula indirectly via driv- 
ing (after /3-decay) the ISM, and stretching the ambient mag- 
netic field, along the jet direction ( Atov an and Dermerl 2003). 
A faster diffusion of the bulk of protons in this direction could 
result in A(E) ^> 0. 1 . Remarkably, in both cases the aspect ra- 
tio should increase at higher energies. 

The HESS collaboration determined the source spectrum 
with the photons from the bright central part of the source 
(within 13.4' from the centre). While our model explains this 
spectrum, it predicts that the energy spectrum will be sig- 
nificantly harder when determined for the entire source (see 
Figure 2). Another specific signature of a GRBR is the un- 
usually hard energy spectrum at E < £ 7l brk ~ 10-100GeV. 
GLAST will be able to verify or falsify this prediction. If 
HESS J1303-631 is a SNR, GLAST will easily detect the 
source at the flux level close to the upper limit of EGRET 
(Mukheriee & Halpern 2005). However, a weak detection 
(confirming the hard spectrum) or non-detection by GLAST 
would signify the tell-tale detection of the signatures of proton 
acceleration by relativistic shocks of a GRB. 

Confirmation of HESS J 1303-631 and at least one other ex- 
tended unidentifi ed TeV source, such as TEV J2032+4130 
JAharonian et alJl2002t lLang et alJl2004l) . as GRBRs would 
confirm the high rate, of order 10~ 4 yr _I , of Galactic GRBs, 
and would make it more likely that GRBs had a direct impact 
on life on Earth ( Melott et alJ2004l) . 
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Space Sciences of the Washington University. 
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